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Abstract 

We study a possible explanation of a 3.0 a excess recently reported by the ATLAS Collaboration in events with X-peaked 
same-flavour opposite-sign lepton pair, jets and large missing transverse momentum in the context of gauge-mediated 
SUSY breaking with more than one hidden sector, the so-called goldstini scenario. In a certain parameter space, the 
gluino two-body decay chain g gxi 2 gZG' becomes dominant, where X 1.2 d" are the Higgsino-like neutralino 
and the massive pseudo-goldstino, respectively, and gluino pair production can contribute to the signal. We find that 
a mass spectrum such as rrig ^ 1000 GeV, m^o ^ ~ 800 GeV and ttiq, ~ 600 GeV demonstrates the rate and the 
distributions of the excess, without conflicting with the stringent constraints from jets plus missing energy analyses and 
with the CMS constraint on the identical final state. 


1. Introduction 

While the LHC Run-I achieved a great success with 
the discovery of the Higgs boson, most of the attempts to 
find new physics failed, and hence we are eagerly waiting 
for results from the LHC Run-H, which has just started 
and entered a new energy frontier. 

However, although no significant signs for new physics 
at the 7 and 8 TeV run have been reported so far, we 
should still keep our eye on the possible faintest imprint. 
The ATLAS Collaboration has recently reported a search 
for supersymmetry (SUSY) with the final state containing 
a pair of same-flavour opposite-sign (SFOS) leptons, jets 
and large missing transverse momentum (^-t) ^ centre- 

of-mass energy of 8 TeV [T]. With data for an integrated 
luminosity of 20.3 fb“^, an intriguing excess of 29 lepton 
pairs peaked at the invariant mass of the Z boson (“on-Z”) 
is observed, while 10.6 ± 3.2 pairs are expected from the 
Standard Model (SM) prediction. The excess corresponds 
to a significance of 3.0 a. 

Interpretation of the excess with some models is not 
a straightforward task, as null results from other searches 
are placing considerably stringent constraints on the vi¬ 
able model parameter space. The ATLAS on-Z signal re¬ 
gion is motivated by search for a pair production of gluino 
in gauge mediation scenarios in the minimal supersym¬ 
metric standard model (MSSM), where gluino decays into 
a neutralino, which subsequently decays into a Z boson 
accompanied by a very light gravitino. The gravitino es¬ 
capes detection, leading to missing momentum. While rel¬ 
atively light gluinos with rrig <1.2 TeV are required [5], 
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studies presented in HE] have shown that explaining the 
excess within this scenario is difficult. At large gluino- 
neutralino mass splitting, gluino decays with a significant 
branching ratio into a top quark, leading to appreciably 
strong bounds from the stop searches. Even in the com¬ 
pressed gluino-neutralino mass region, the process with 
the hadronic decay mode of the energetic Z boson is con¬ 
strained by the dedicated jets plus $t searches. Given 
that situation, several solutions have been proposed in 
various models, e.g. NMSSM [31 E], MSSM with a light 
sbottom or stop |5], pMSSM |7], split SUSY j^, as well as 
non-SUSY models liiini. 

In this paper, we investigate an alternative SUSY model, 
i.e. a model with multiple hidden sectors, the so-called 
goldstini model [TTUI4] . especially in the context of gauge 
mediated SUSY breaking. We consider a gluino (g), the 
lightest and second lightest neutralinos (xi_ 2 )j a pseudo- 
goldstino {&) in the spectrum as a simplified model. As 
we verify below, in a certain parameter space, the gluino 
two-body decay chain 

9^9 + Xi,2 9 + Z + G' (1) 

becomes dominant, and gluino pair production can con¬ 
tribute to the signal. Jets from the gluino decay and the 
hadronic Z decay are softened when the mass spectrum 
is compressed due to the massive nature of the pseudo- 
goldstino. We find that there is a viable parameter space 
even after the multijet plus constraint m. as well 
as the GMS constraint on the identical final state [16] . 
are taken into account. We also show that the two-body 
gluino decay in 0 provides a better fit to the data for the 
distributions with respect to the three-body gluino decay 
9 qqxi,2- 
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The paper is organized as follows: In Sec. we briefly 
review the model involving two SUSY breaking sectors, 
and verify the gluino two-body decay chain in Q. In 
Sec. we recast the ATLAS on-Z analysis for our bench¬ 
mark scenarios to find a viable parameter space. The 
constraints coming from other searches are also discussed. 
Sec. 1^ is devoted to our summary. In [Appendix A| we 
give the explicit formulas for the neutralino decays, while 
Appendix Bjwe present the validation of the implemen¬ 


tation for the analyses. 
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2. Model 

The model we consider, in order to explain the AT¬ 
LAS Z-peaked excess, is a scenario of gauge mediation 
with more than one secluded SUSY breaking sector. Each 
SUSY breaking sector contains a massless goldstino, i.e. the 
goldstone fermion of spontaneous SUSY breaking. In a 
two-sector scenario there are then two massless goldstini. 
Models of goldstini have been studied in [nHH iizHini 
HOHin], and the description of goldstini in gauge media¬ 
tion models and the quantum corrections to the pseudo- 
goldstino mass have been investigated in [14]. Indeed, 
when taking into account the mediation mechanism and 
the interaction of the SUSY breaking sector with the MSSM, 
a linear combination of the goldstini becomes the true 
goldstino (the longitudinal component of the gravitino), 
while the orthogonal combination is a pseudo-goldstino. 
Denoting with fa and Ga the SUSY breaking scales and 
the goldstini of the two sectors respectively, the true (mass¬ 
less) goldstino and the pseudo-goldstino are 

G= |(/iGi+/2G2), G'= |(-/2Gi+/iG2), (2) 

where F = \Jf\F f\ is the total SUSY breaking scale. 
The true goldstino (or gravitino) mass is related to the 
SUSY breaking scale and the Planck scale as vHq cx F/Mpi, 
and hence the low-energy SUSY breaking scenario as in 
gauge mediation leads to a very light goldstino. On the 
other hand, the mass of the pseudo-goldstino is not pro¬ 
tected by any symmetry and generically receives relevant 
quantum corrections, proportional to the SUSY breaking 
terms El- In our paper, we consider the pseudo-goldstino 
mass to be of the order of a few hundreds of GeV. 

The interaction of each of the goldstini with the par¬ 
ticles in the visible sector are determined by the super¬ 
current, and hence the couplings of the true goldstino and 
of the pseudo-goldstino with the MSSM states are also 
fixed in terms of the soft SUSY breaking masses. The 
noteworthy feature is that the pseudo-goldstino couplings 
are generically different with respect to the goldstino cou¬ 
plings. Indeed they depend on the relative contribution to 
the soft masses from each of the SUSY breaking sector, and 
they can be enhanced without changing the overall soft 
masses. This implies that the decay of SUSY particles can 
be dominantly into a massive pseudo-goldstino plus a SM 


Figure 1: Mass spectrum for our simplified goldstini model, with 
the relevant decay modes. 

particle, drastically changing the final state topology with 
respect to the usual decay into the true goldstino [20H25] . 

In this paper we consider a simplified model including 
as low energy degrees of freedom only a gluino, Higgsino- 
like neutralinos, a pseudo-goldstino and a goldstino. The 
mass spectrum of the model is depicted in Fig.j^ Note that 
the Higgsino fields include two almost degenerate neutral 
mass eigenstates and a charged one. We also note that the 
true goldstino is in the bottom of the spectrum, but is not 
shown in Fig. since it is irrelevant in our scenario. In the 
following we study each decay step in detail to verify the 
gluino decay chain in Q. 

Gluino decay. In the spectrum presented in Fig. the 
gluino can potentially have several decay channels, de¬ 
pending on the mass splitting between the gluino and the 
neutralinos. 

In case of the large splitting such as {rrig — m^o) ^ 
mt, the tree-level three-body decays into a pair of third- 
generation quarks and a chargino/neutralino are domi¬ 
nant |4]. 

On the other hand, in the regime 

Amg_^o ^ = rrig — m^o ^ < 200 GeV, (3) 

the gluino decays predominantly into a gluon and a neu¬ 
tralino via a (s)top loop. The analytic expression for the 
9 9Xi decay can be found in PSI E7] . We find that this 
result is robust as soon as the Bino (Mb) and Wino (Mw) 
masses are moderately larger than the Higgsino mass (p). 
We checked this decay pattern with SUSY-HIT |55|, for 
instance, by fixing nig = 1000 GeV, fj. = 800 GeV, and 
tan/3 = 10, we find B{g —)• gxi,2) > 85% as soon as 
Mb and Mw are larger than about 1.4 TeV, with squark 
masses of the order G(5) TeV. 

The gluino decays into a pseudo-goldstino or a gold¬ 
stino with a gluon is also possible. However, unless the 
gluino coupling to the pseudo-goldstino is largely enhanced, 
these decay modes are always suppressed compared with 
the decays into the MSSM states. 

In this paper, to fit the ATLAS excess, we consider the 
small mass splitting in Eq. where the only g —> gffl 2 
decay is relevant. 
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Neutralino decay. The interaction lagrangian which is rel¬ 
evant for the neutralino decay into a pseudo-goldstino and 
a electroweak gauge boson or a Higgs boson is given by El 

naiiiHis] 






2y/2F 


. Vz-r 
''2^/2F 


G'a^a^'x^Z^ 


Vk 

y/2F 


X°G'h + h.c., (4) 


where and are the field strengths of the photon 
and the Z boson, respectively, and h is the lightest Higgs 
boson in the decoupling limit. The goldstino lagrangian 
is the same, but with different coefficients, y ^ y. As 
mentioned above, the pseudo-goldstino couplings y can be 
larger than the goldstino couplings y, and in a simplified 
model approach they can be considered as free parame¬ 
ters [ 2 TI [ 22 I [ 2 I] . Given a set of soft terms originating from 
the two SUSY breaking sectors, one can compute these 
couplings with the formulas collected in [Appendix A[ 
Although the neutralino decay can present a rich pat¬ 
tern, with six competing decay modes as {Z, h, 7 ) plus 
(G',G), we are interested in a scenario where the neu¬ 
tralino predominantly decays into a pseudo-goldstino and 
a Z boson. This scenario can be realized by enhanc¬ 
ing the coupling parameters y, especially yz^ and yzL^ 
and by assuming the Higgsino-like neutralino. The de¬ 
cay formulas are reported in [Appendix A where we also 
provide the details of our illustrative benchmark point in 
the SUSY breaking parameters determining the couplings 
y and y; we typically take y ~ 800 GeV and Mb = 
Mw ~ 1-5 TeV. For the parameters we consider in this 
paper, the mass splitting between the two neutralinos is 
of the order of a few GeV. In this scenario, the second 
lightest neutralino y® decays to the lightest neutralino 
Xi with soft SM particle emissions. Indeed we checked 
with SUSY-HIT [55] and the formula in Appendix A 


that its decay modes to the goldstino and to the pseudo- 
goldstino for our benchmark point are negligible compared 
with the X 2 Xi decays. Hence in the following we as¬ 
sume B{x 2 Xi + undetectable SM particles) = 100%. 
The final state topology is then determined by the possi¬ 
ble Xi decays, which we now investigate. 

In Fig. we show the branching ratios of the light¬ 
est neutralino as a function of the pseudo-goldstino mass, 
evaluated on our illustrative benchmark point described in 
detail in [Appendix A| The decay pattern depends on the 
mass splitting Am^o_Q, = m^o — m^, and on the possi¬ 
ble kinematically allowed modes. The branching ratio of 
Xi —>■ ZG' is always greater than 80% for Am^o_Q, > ruh, 
and saturates at 100% for niz < Am^o_Q, < ruh- We note 

that the Xi iG' decay is negligible in our parameter 
choice. In the regime of Am^o_Q, < mz, the decays into 
a true goldstino plus a 7 , Z or h become dominant due to 
the phase space. In the following, we consider the region 

Am-o_( 5 , > (5) 



Figure 2: Branching ratios of the lightest neutralino as a function 
of the p seudo-goldstino mass for a representative benchmark point 
given in 
the plot. 


Appendix A The Xi 1 ^' decay is too small to see in 


where the decay is dominant, as assumed in the 

simplified model in Fig. and shown in Fig. 

We have also verified that the total decay width in 
the region of interest is always larger than 2 x 10“^^ GeV, 
implying that the decays happen promptly in the detector. 

Finally, the pseudo-goldstino will eventually decay into 
a goldstino plus 7 , Z or h. However, one can verify that 
the pseudo-goldstino is enough long-lived so that the decay 
happens outside the detector. With the formulas listed 
in [Appendix A[ one can indeed compute the pseudo- 
goldstino decay width and we find that for the benchmark 
point it is around 10“^^ — 10“^^ GeV, i.e. Tq, < 1 sec, 
depending on the pseudo-goldstino mass. Even though we 
are not addressing cosmological issues in this paper, we 
observe that this decay is fast enough not to spoil the big 
bang nucleosynthesis (BBN) [55] . 


3. Analyses 


At the LHG, our simplified goldstini model can be 
probed by gluino production. Given the discussion on the 
decays in the previous section, the gluino pair production 
leads to the process illustrated in Fig. [^ where we assume 
the 100 % branching ratio at each decay step as a good 
approximation. Depending on the decay of the Z boson, 
the final state can be SFOS lepton pair -|- jets -I- and 
jets 4- 

In this section, we investigate whether or not the event 
topology of our model can successfully fit the ATLAS on- 
Z excess without conflicting with other searches. The most 
stringent constraints come from ATLAS multijet search m 
and CMS dilepton search [T5] in the parameter region of 
our interest [31115]. The latter has signal regions which 
look at the same final state as in the ATLAS on-Z signal 
region and potentially quite constraining. We have imple¬ 
mented these analyses as well as the ATLAS on-Z signal 
region in Atom |30| . Some description and validation re¬ 
sults of Atom are given in Appendix B 
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Figure 3: The process at the LHC in our simplified goldstini model. 


In order to fit the excess we scan the gluino mass having 
the neutralino masses fixed at 


TO-o ^=m~g- 200 GeV. (6) 


We consider three cases featuring the pseudo-goldstino 
mass: 


rriQ, 


0 (A) 

m^o - 200 GeV (B) 
m^o - 100 GeV (C) 


(7) 


Case A is equivalent to the gauge mediation scenario with 
only one SUSY breaking sector and a very light gravitino, 
while cases B and C have compressed spectra. 

In order to assess the consistency between the model 
and data, the ideal approach would be to carry out a global 
fit, treating the excess and constraints in the same manner. 
This requires the details of the systematic uncertainties 
and good understanding of the correlation among different 
signal regions. Rather than taking this rigorous approach, 
in this exploratory paper we instead fit the model to the 
excess independently from the constraints and check the 
exclusion individually for the signal regions using the fol¬ 
lowing prescription. 

To see a goodness of the fit, we define the measure R 
as 

R = Nsvsy/ (A^obs — A^sm), (8 ) 

for the ATLAS on-Z signal region, where Vsusy is the 
expected SUSY events, Vobs is the number of observed 
events and Nqm is the expected SM events in the signal 
region. With this definition the best fit is given at R = 1. 

For the other signal regions, labeled i, used as con¬ 
straint, we instead define 

R ' ^ AT^usy/A^bsm. (9) 


where Vsusy is ttie expected SUSY events and A^egM i® 
95% CLs limit obtained from signal region i. Having any 
Ri greater than one indicates that the model is strongly 
disfavoured. 


The VgusY can b® expressed as agg ■ C ■ e(q . Here C is 
the integrated luminosity used in the analysis and agg is 
the production cross section of the gluino pair, for which 
we use the values reported in [SH |32j . To estimate the effi¬ 
ciency e(i) we use the following simulation chain: first the 
signal events are generated using MadGraph5_aMC(§NL0 [33] 
and showered and hadronized by PythiaS |34] . The hadron 
level events are then passed to Atom m to estimate the 
efficiency for each signal region taking the detector effects 
into account. For our signal simulation, we extended the 
goldstini model |3T1 |3Sj (building on [53] ) to include the 
two-body gluino decay by using FeynRules2 m- 

The results of the three cases of the goldstini scenario 
in Eq. Q are presented in Figs. |^and[5j where the modes 
are confronted with the ATLAS multijet [T5| and the CMS 
dilepton [T3] searches, respectively. The fitting measure R 
for the ATLAS on-Z signal region [T] is shown with the 
solid black curve, whereas the constraints i?* are shown 
with the other curves, corresponding to the signal regions 
(2jl, 2jm, 2jt, 3j, 4jl-, 4jl, 4jm, 4jt, 5j, 6 j, 6 jm, 6 jt, 6 j-P) in 
the ATLAS multijet search in Fig. I^and (cms 2 jl, cms 2 jm, 
cms2jh, cms3jl, cms3jm, cms3jh, cms(G), cms(F)) in the 
CMS dilepton search in Fig. jS^I The green (yellow) band 
around i? = I represents the I (2) a region of the fit for 
the ATLAS on-Z excess. 

One can see that the entire region of case A, i.e. the 
massless goldstino case, is excluded by the ATLAS multijet 
search. The CMS dilepton search also excludes the region 
with rrig < l.I TeV. These strong limits can be attributed 
to the large mass gap between the gluino and the pseudo- 
goldstino, because of which the jets and leptons from the 
gluino decays tend to be hard, helping to increase the ef¬ 
ficiency of the ATLAS multijet search as well as the CMS 
dilepton search. 

Moving to case B and C, one can see that the con¬ 
straints are generally more relaxed than in case A, since 
these cases have milder mass hierarchy with the massive 
pseudo-goldstino. The ATLAS multijet search now ex¬ 
cludes the gluino mass up to 970 (810) GeV for case B 
(C). Here, the best fit of the ATLAS on-Z excess is given 
at rrig = 980 (900) GeV for case B (C), which is just out¬ 
side of the ATLAS multijet exclusion limit. 

The CMS dilepton search, on the other hand, still pro¬ 
vides tight constraints, especially the cms 2 jh signal re¬ 
gion, where > 300 GeV is required, very similar to 
the > 225 GeV cut in the ATLAS on-Z analysis. 
The difference between the two analyses is that the AT¬ 
LAS search additionally imposes Ht > 600 GeV, where 
Ht is the scalar sum of the px of jets and leptons. For 
case B, the best fit point {rrig = 980 GeV) is just on the 


^ The nj in the signal region name indicates it requires more than 
n high pt jets. The letter “1”, “m”, “t” for the ATLAS multijet 
search means “loose”, “medium”, “tight”, while “1”, “m”, “h” for 
the CMS dilepton search denotes “low”, “medium”, “high”. The 
cms(C) and cms(F) represent the central and forward signal regions 
in which the number of jets and are treated inclusively. See m 
and m for the exact definition. 
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Figure 4: R values for the fit (black solid) and constraints (others) for cases A (left), B (middle) and C (right) in Eq. 0- The green and 
yellow bands correspond to the 1 and 2 a regions of the fit. If the black curve is in the bands, the model provides a good fit for the ATLAS 
on-Z excess On the other hand, if there is any other curve above one, the model point is strongly disfavoured by the signal region 
corresponding to the curve in the ATLAS multijet + analyses [15| . 



Figure 5: The same as in Fig.j^ but for the constraints from the CMS on-Z analyses m- 
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Figure 6: Comparison between data and signal + background in the (left) and jet multiplicity (right) distributions in the ATLAS on-Z 
signal region [J at our best fit point: ( 77 X 3 , 771-0 ,m^,) = (1000,800,600) GeV. 
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Figure 7: Comparison between data and signal + background in the distribution in the CMS on-Z signal region with > 2 jets (left) and 
> 3 jets (right) 1161 at our best fit point: (mg,m-o ,m^,) = (1000,800,600) GeV. 


exclusion boundary, while for case C the best ht point 
{rrig = 900 GeV) is just excluded by the cms2jh signal re¬ 
gion. For case B, at rrig = 1 TeV, the fit for the ATLAS 
on-Z excess is still within the 1 a band, and the point is 
still outside of the 95% CLg exclusion. We therefore choose 
our best ht benchmark point as 


(rug, TO^o m^,) = (1000,800,600) GeV (10) 

for the following analysis. Even for case G, the tension be¬ 
tween the data and the prediction observed in the ATLAS 
on-Z signal region can be ameliorated to the 2 cr level with 
rrig = 950 GeV, which is outside the 95% GL^ exclusion 
region. 

In Fig. we compare the data with the signal -|- back¬ 
ground in the (left) and jet multiplicity (right) dis¬ 
tributions in the ATLAS on-Z signal region at our best 


ht point in (10). Here we took the data and the SM 


background from Figs. 6 and 7 in the ATLAS paper [I] 
and combined ee and channels. The data in the 
distribution has a preference of low peaking around 
240 GeV and roughly falling down up to around 500 GeV. 
The distribution is well htted with the signal -I- back¬ 
ground at our best ht point because the massive nature 
of the pseudo-goldstino reduces the in the event. In 
the jet multiplicity distribution the data prefers 2 — 5 jets 
and disfavors the region with > 6 jets. The distribution of 
the signal -I- background at our best ht point peaks around 
3 — 4 and gives a good ht to the data. This is an advan¬ 
tage of the radiative decay g —^ gx)* 2 compared to the 
three-body g —>■ qqxi 2 decay because the number of jets 
is reduced typically by two. 

As a reference, in Fig.j^we also compare the signal with 
the CMS data in the on-Z signal region with Ajets > 2 and 
> 3 at our best ht point in ( [To| ) . Here we took the data and 
the SM background from Fig. 7 in the CMS paper [16]. As 
already seen in Fig.[^ the most stringent constraint comes 
from the high region, i.e. cms2jh and cms3jh, where 
the $x > 300 GeV bin is considered. 


4. Summary 

Even though LHC at 8 TeV has not discovered a new 
physics signal, the Run-I data contain a few small excesses 
that deserve a thorough investigation. Recently the AT¬ 
LAS Collaboration has reported a 3.0 a excess in dilepton 
-I- jets -I- dilepton reconstructing the Z-boson 

mass. 

In this paper we have proposed an explanation of such 
excess in a SUSY model of gauge mediation with two 
SUSY breaking sectors, presenting in the SUSY spectrum 
an extra neutral fermion besides the MSSM neutralinos, 
that is the pseudo-goldstino. Our simplified model con¬ 
sists in a gluino, a pair of Higgsino-like neutralinos, and 
a pseudo-goldstino. We showed that our goldstini model 
can explain the ATLAS .Z-peaked excess without conflict¬ 
ing with the constraints from multijet -I- "well as from 

the CMS analysis for the same final state. The compressed 
spectrum such as rrig ~ 1000 GeV, m^o ^ 800 GeV and 
rUg, ~ 600 GeV gives a very good fit to the data not only 
for the rate but also for the distributions. 

If this is indeed the origin of the excess, the 13-TeV 
LHC with both the ATLAS and CMS Collaborations will 
soon confirm the excess with much greater significance. 
We await the verdict from the LHC Run-H. 
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Appendix A. Neutralino decay 

In this appendix, we review the true goldstino and the 
pseudo-goldstino couplings with the MSSM states, and de¬ 
rive the formulas for the relevant decay widths; see also 
Refs. [II1II91EIH23]. 

In the two-sector goldstini model, the lagrangian rele¬ 
vant to the neutralino decay, in the gauge eigenstate (i/'i = 

CD - - PaiGaX^ - \j^abGaGb - 

+ TaiGa-iptho + GijPpla^^PpjZ^ + iLia-ipia>^a''GaZ^^, 

(A.l) 


where Mf, is the standard 4x4 neutralino mass matrix. 

•'J 

The other parameters are 


Pa = 




Ta = 


V2f, 

( 


V 




9lSfi 

-92313 

0 

0 


MB{a){DY) 
Mw(a){DT3) 

-H, 

V V2v{mj^^^^^Si3 - B^a)Cl3) ) 

0 0 -gicp 

0 0 52C/3 

-fflC/3 92Cp 0 

fflS/3 -525/3 0 0 / 

( mzMB{a)SwC2p \ 


G = 


L„. = 


V^fa ''^Hd(a)'^P ^{a)Sp 

V ~ B(^a)Cp 

/ 0 0 0 0 \ 

92 0 0 0 0 

2 cw 0 0 10 

\ 0 0 0 -1 y 

/ —MB(a)SW ^ 
Mw{a)CW 

2V2fa 0 

V 0 


I 


(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 


where we denote with MB/w(a)i (a) 

contribution to the soft term from sector a, and consider 
the decoupling limit for the Higgses. For convenience it is 
also useful to define tilded soft masses, e.g. 


Mb 


h 

fi 


£ 

Mb{1) + -wMb(2), 
J2 


(A.7) 


while the non-tilded are the physical ones, e.g. 


Mb = Mb{i) + Mb{2)- 


The mass matrix is 


M^ 


-f,Ml2 Mi2 \ 

Mi2 -%Mi2 J ’ 


(A.8) 


(A.9) 


and it is such that it has one zero eigenvalues along the 
true goldstino eigenstate in Eq. ([^. The entry AI 12 is a 
model dependent quantity that has been computed in [14] . 
The other eigenvalue of M^ is the pseudo-goldstino mass 
m-G" 


Given a set of soft terms and the relative contributions 
from the two sectors, one should diagonalize the lagrangian 
and write the resulting couplings in the mass eigenbasis. 
This ends up in the lagrangian quoted in Eq. Q, whose 
couplings can be computed for a given set of soft terms. 
Note that, even though suppressed by l/fa, there will be 
some mixing between the neutralinos and the goldstini. 
This is the reason why we have also written in (A.l) the 
terms containing only the neutralino gauge eigenstates. 
For instance, once rotated into the mass eigenbasis, these 
interactions generate the effective couplings yz^ in Eq. (|^. 

Although we study a rather heavy pseudo-goldstino 
case in this paper, we consider the rriQ, = 0 limit for 
a while to get an insight on the feature of the pseudo- 
goldstino couplings by the analytic expressions. At leading 
order in razlra^ in the neutralino and pseudo-goldstino 
mixing and neglecting the pseudo-goldstino mass, the ef¬ 
fective couplings in the pseudo-goldstino lagrangian 0 are 


= m.^o{KBNiiCw + Kw^asw) 

+ mz{N*^sp — N*^cp)swcw{Kb — Kw), (A.10) 
Vzt — m^oi—KBNiisw + KwNi2Cw) 

+ Tnz{N*j^sp — N*^cp){s‘^Kb + c^Ky/), (A.11) 
yZi = —TnzK^{—NiiSw + Ni2Cw) 

-m^o{KuN*^sp-KdN*^cp), (A.12) 

Vh = —wz cos2/3{—KbMbN*iSw + KwMwN*2Cw) 
-p‘^{KdN*3Cp + KuN*^sp), (A.13) 


where the neutralino rotation matrix Nij is defined as 
in EHl- The K factors read 


Kb = 


Mb 
Mb' 
1 


Kw = 


M, 


w 


Mw ’ 


+ 3^Ku 


- Stan/ 


H—^ {sw^B + CwKw ) cos 2/3), 

Ky, = - \ (rfiB - Boot 13 
2 

-^iswKB + c^Kw)cos2l3), (A.14) 


^We follow the conventions of m- 
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where we keep the terms of order in the last two K 
factors in order to show the goldstino limit explicitly. The 
goldstino lagrangian is recovered by converting all 

the tilded soft terms to the un-tilded ones, and using the 
electroweak symmetry breaking (EWSB) conditions 


0 = +/i^ — Bcot/3 —cos2 /?, (A. 15) 

0 = — Stan/? + cos2/3, (A. 16) 

2 


resulting in all the K factors in Eq. (A.141 to be unity. 
From the above analytic couplings, for instance, one can 
observe that for a mostly Higgsino-like neutralino with 
Kb = Kw the pseudo-goldstino coupling to the photon 
is suppressed. We also checked that the above expressions 
agree with the numerics if the pseudo-goldstino mass is 
sufficiently smaller than the neutralino masses. For the 
large pseudo-goldstino mass, on the other hand, there can 
be deviations from these formulas once we rotate in the 
mass eigenbasis. Hence, in order to compute correctly the 
effective couplings in Eq. Q for generic pseudo-goldstino 
masses, one is instructed to diagonalize numerically the 
lagrangian (A.l), for a given set of soft terms, and write 
it in the mass eigenbasis. 

Once the effective couplings in Eq. Q are evaluated, 
the neutralino decay widths can be computed as 


^X°^ZG' 


^hG' 


IGttF^ TO-o/ 


ml' 

32nF‘^ \ m^o) rniq 
X +mG,f+ ml) 


+ + ^G'f + 2 w|) 

+ QyzTyz^mzim^o + TOq,)] , 

Ph{yl?m^o / \ 

327rF2 m| J ’ 


(A. 17) 


(A. 18) 
(A. 19) 


where ^z/h = /3(TOg,/w|o, with,5(a, 6) = (1 + 

+ 6^ — 2a — 26 — 2a6)^/^. The standard decays into the 
massless goldstino are obtained from this expression by 
sending the tilded quantities to the non-tilded quantities 
and putting TOq, —>■ 0 . 

In this paper we are interested in configurations where 
the two lightest neutralinos are mostly Higgsinos and the 
decay of is predominantly into a massive pseudo-goldstino 
and a Z boson. To find if this scenario is compatible with 
the lagrangian just explained, we numerically explored the 
parameters of the model (i.e. the soft terms) looking for 
a representative benchmark point satisfying these require¬ 
ments. In Fig. we report the branching ratios for a con¬ 
figuration with fj, — 804 GeV, \/B = 800 GeV, Mb = 
Mw = 1.5 TeV and tan/3 = 10. The total soft masses for 
the Higgses are extracted by solving the EWSB conditions. 


The two SUSY breaking scales are chosen as v7i = 1-5 x 
10® GeV and v72 = 5 x 10^ GeV. The gaugino masses and 
the Higgs soft masses are distributed in the two sectors as 
Mb/w{i)IMbiw{ 2 ) = tan‘^0, = cot^ff. 

The angle 0 is taken to be tan0 = 0.2. Finally the B 
terms are chosen as B(^a) = fa/ifi + f 2 )B. With these val¬ 
ues we obtain a scenario where the two lightest neutralinos 
are mostly Higgsinos, with masses at (797,805) GeV. The 
lightest neutralino decay is predominantly into a pseudo- 
goldstino plus a Z boson for different ranges of pseudo- 
goldstino mass, as we show in Fig. The pseudo-goldstino 
mass is varied by changing AI 12 consistently with the per¬ 
turbative definition in M- The decay width of the neu¬ 
tralinos cited in the text are computed with the decay 
formulas (A. 171, (A. 18) and (A.19). Moreover, the de¬ 
cay of the pseudo-goldstino can be computed using the 
same formula quoted above, where the coupling between 
the pseudo-goldstino and the goldstino are extracted nu¬ 
merically from the original lagrangian once we switch to 
the mass eigenbasis. 


Appendix B. Analysis implementation 

For this paper, we implemented several ATLAS and 
GMS analyses in the Atom framework [30]. Atom is a pro¬ 
gram based on Rivet [43] and maps the truth level par¬ 
ticles into the reconstructed objects such as isolated elec¬ 
trons and 6-jets according to the detector performances 
reported by ATLAS and GMS. The validation and some 
application of the code can be found in |44]-[48j. 

In Table El] we show some of the validation results 
as an example. The numbers in the second column rep¬ 
resent the expected signal events for each step of the cut 
used in the 5j signal region (SR) reported by the ATLAS 
multijet analysis [15], based on the ql —>■ qxf 
topology with (m,, m-±, Wj^o) = (665,465,265) GeV. The 
right column shows the ratios between the Atom and AT¬ 
LAS results. One can see that these ratios are close to 
one within ^ 20% accuracy indicating a good agreement 
between the Atom and ATLAS simulations. 


5j SR cuts 

-''SUSY 

Atom/Exp 

> 160, Pr'"’ > 130(60) GeV 

317.3 

1.17 

> 60 GeV 

306.2 

1.12 

> 60 GeV 

247.6 

1.04 

p^ > 60 GeV 

141.8 

1.00 

^</’(jl,2,3 ) ^t) > 0.4 

118.6 

1.01 

A</.(A >3 >40GeV,.^T) >0.2 

103.1 

1.01 

$T/'mes{Nj) > 0.2 

85.6 

1.04 

TOeff(incl.) > 1200 GeV 

20.5 

1.18 


Table B.l: “5j” SR validation table for the implementation of the 
ATLAS multijet analysis in Atom. The decay chain for validation in 
consideration is —>■ qx^ —>■ qW^x^. 
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